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Abstract. Granular segregation is ubiquitous in industrial, geological or daily-life 
context, but there is still no unifying theoretical approach. In this review, we examine 
two examples of granular segregation - shallow rapid flows and rotating drums - 
which suggest that the dynamics of systems at intermediate and high density might 
be amenable to a statistical mechanics approach. 
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He wondered how long it would take before everything in the bowl was completely 
mixed. But what did that mean, completely mixed? Every grain of the 
ingredients would have to be distributed perfectly, the particles of salt and baking 
soda spaced just so throughout the flour, each one a fixed distance from all the 
rest. He tried to picture it, a solid, three-dimensional white field supporting a 
dense and uniform lattice of particles of other shades of white. And what about 
the flour itself, no two grains of which were alike — how could that be completely 
mixed, even if there were no other ingredients present in it, making their own 
pattern? And how would he know when that moment of perfect distribution had 
been achieved? 

John Banville, The Infinities 

As even fictional characters have come to notice, the mixing of granular materials 
is a tricky business. The very processes by which you might hope to evenly distribute 
particles of different type - stirring or shaking - quite commonly, in fact, create domains 
of different particle types [U [2]. The aim of this review is to explore the problem 
of mixing and segregation in dense granular systems within the context of statistical 
mechanics, and then to use this approach to highlight open research questions. 

We begin by reviewing why ordinary fluids either mix or phase-separate, and then 
discuss the extent to which granular materials exhibit similar and/or different behaviors. 
A familiar example is the success of kinetic theory in describing Soret-like effects in 
granular gases. We follow with two examples from higher-density systems which have 
been more challenging to understand. First, we examine segregation dynamics in shallow 
flows, where recent theories about the local distribution of particles and voids [31 H] 
provide hope for new approaches. Second, we suggest that both axial segregation and its 
subsequent coarsening observed in rotating drums can be understood as a surface tension 
driven instability within the statistical mechanics framework suggested by Edwards and 
coworkers [6] . 

1. Segregation: ordinary vs. granular fluids 

Classical statistical mechanics - in the form of thermal fluctuations, entropy, free energy, 
and gradients - has allowed us to understand why and how a mixture of oil and 
water will spontaneously phase-separate, while a mixture of vinegar and water will 
not. In equilibrium thermal systems (at fixed temperature and volume), the final state 
is governed by the minimization of free energy T = U — TS, and thermal fluctuations 
allow the system to explore many different configurations until the minimum is reached. 
For conventional fluids, the molecular sizes are typically similar enough that depletion 
effects are small; in such cases, the entropy 5* will favor a mixed state. Therefore, the 
presence or absence of phase separation will largely be a consequence of minimizing the 
potential energy U . If we consider two species 1 and 2, the three contributions to U are 
the interaction energies Uu, U 2 2, and Lq 2 - If the interaction between unlike molecules, 
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U12, is less negative (smaller magnitude) than U\\ and U22, then the potential energy 
contribution will favor phase separation. If AU = U\ 2 — (Uu + U^)j1 is large enough 
to overcome the mixing entropy, the equilibrium state is a segregated one. 

Even if the system should, by these arguments, de-mix, the surface tension a 
between phases 1 and 2 can nonetheless cause the segregation process to be kinetically 
hindered. A droplet of size R containing a single species loses free energy via the creation 
of an interface of surface area A, by an amount AJ 7 = —Aa oc R 2 . Simultaneously, it 
gains free energy proportional to its volume V, in an amount AJ 7 = V AU oc R 3 . A 
nucleated droplet above a certain critical size R c will have a net negative change in free 
energy, and will therefore continue to grow in size; droplets smaller than R c will shrink. 
Therefore, a phase-separated state can only occur when thermal fluctuations provide at 
least one nucleated droplet above this critical size. In the absence of such a droplet, the 
system can stay in the metastable mixed state indefinitely. 

When a thermal system is driven out of equilibrium, whether by gravitational, 
thermal, or other gradients, differences in the transport coefficients can provide 
additional sources of demixing. A notable example is the effect of thermal gradients 
on mixtures of ordinary fluids. A phenomenon known as thermophoresis [8] breaks 
the symmetry such that molecules of differing masses migrate either up or down the 
thermal gradient. Typically, the heavier species movies towards the lower-temperature 
region, commonly known as the Soret effect. While thermophoresis is exploited 




Figure 1. Thermal diffusion in a granular gas. (a) Image of a vertically shaken, 
maximally expanded mixture of equal volumes of brass spheres with 1.19 mm and 
2.38 mm diameter. Shaking parameters are 15 Hz and an average acceleration of 7g. 
The larger particles accumulate in a small region in the lower half of the column, 
(b) Results from a MD simulation which matches the experiments shown in panel a. 
The granular temperature shows a pronounced minimum in which the large particles 
preferentially accumulate. This agrees with the prediction of granular kinetic theory 
that the larger particles have a higher thermal diffusion coefficient. Modified from [7]. 
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in industrial processes, and thought to be an important factor in transport across 
biological membranes, the mechanisms underlying this behavior remain a matter of 
some debate [9]. In practice, the Soret coefficient remains an empirically-measured 
transport coefficient. 

It is not a priori clear how many of these explanations will survive in the case of 
granular systems, which are formed out of macroscopic constituents. First, the large 
size of the particles means that the thermal energy k^T is typically a factor of 10 12 
smaller than other energy scales. As a consequence, granular systems are said to 
be athermal, and any dynamics are the consequence of external driving rather than 
thermal fluctuations. Second, due to both friction and inelastic collisions, the system is 
dissipative: a steady energy influx is required for any persistent dynamics, and hence 
any mixing or segregation of particles. 

In spite of these caveats, there is hope that a statistical mechanics-like approach 
can elucidate granular dynamics. In the case of dilute granular gases, segregation can be 
understood on the basis of an appropriately-extended kinetic theory. This approach first 
defines a granular temperature T g as the averaged kinetic energy of the random motion 
of the particles [TO]. Within this framework, it is possible to compute the dependence 
of the transport coefficients on the particle properties. Because most granular systems 
are driven by injecting energy at the boundaries, gradients in T g are ubiquitous. An 
example is shown in Fig. [T], in which the large particles accumulate at a well-defined 
temperature minimum. Such Soret effects have been well-described by granular kinetic 

theory [in H21 H3i nn nn H71 hei unj. 

2. Vertical segregation in rapid, shallow flows 

There is a common intermediate regime of granular flows, where particles are freely- 
flowing but nonetheless experience significant constraints on their motion due to the 
nearness of neighbors. In such cases, the assumptions of kinetic theory have broken 
down. Such flows are easily found in both natural systems (debris and pyroclastic flows 
[20]) and many industrial processes [2"Tj . 

Such flows are commonly shallow, and large particles segregate to the free upper 
surface of the flow. Therefore, theoretical approaches have paid particular attention to 
this geometry. Models by a number of groups [22j E31 EH EHl e.g.] rely on the idea 
of kinetic sieving as the source of size-segregation. In each case, the underlying cause 
lies in the comfortable and oft-stated idea that as sheared layers flow past each other, 
small particles are more likely to find voids to fall into than are the large particles [26] . 
This imbalance in the downward flux of small particles as compared to large particles, 
combined with a continuity equation for the conservation of mass, will generically lead to 
size-segregation. However, this sieving effect has yet to be derived from first principles. 

Savage and Lun [22] provided an early attempt at a quantitative model, which they 
based on a very general maximum-entropy approach. In this approach, a particle has a 
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Figure 2. Schematic drawing of segregation within an shallow flow. An initial mixture 
with concentration £ = ~ (at all depths) evolves to one in which £ ss at the upper 
surface and £ « 1 at the base of the flow, after exposure to a flow velocity v = (u, v) 
which generates a segregation velocity v s via kinetic sieving. 



probability V(e) of finding a void of lateral size e 

V{e) = ^— exp (1) 

E — Eq |_ £ ~ £ 0_ 

with e the mean void size and eq the minimum void size. This probability is then 
used to calculate the probability of one layer of the flow capturing a particle from 
a neighboring layer. The resulting model provides a prediction for the steady-state 
particle size distribution in a steady shear flow. 

One of the shortcomings of the Savage and Lun| model is that it predicts segregation 
even in the absence of gravity. The model of Gray and Thornton solved this difficulty by 
introducing momentum-balance in deriving an equation for the evolving concentration 
field £(x, y, t) of the small particles, where £ = 1 (0) corresponds to all small (large) 
particles. For the situation shown in Fig. |2j 

where v = (u, v) is the velocity field of the bulk granular material in the x = (x, y) 
directions, as determined from an either experimental measurements or a rheological 
model. A segregation velocity f s (0 is superimposed on the small particles, relative to 
the bulk motion v. In Gray and Thornton [21], the functional form is taken to be 
v a (£) = S(l — £), so that the last term of Eq. [2] is nonlinear. This is a lowest-order 
(linear perturbation) model in which the large particles (which have a concentration 
1—0 experience a pressure-imbalance with respect to the small particles, and S is 
the non-dimensional strength of the segregation velocity. While the model has recently 
seen some semi-quantitative success in predicting the dynamics of both experiments 
and DEM simulations, [271 [28] , the correct functional form of v s (£, <p, P, ■ ■ ■) remains 
little-investigated. 
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Figure 3. (a) Sample image showing detected voids (local e), using previously- 
published dataset [35] of an agitated two-dimensional floating on a horizontal air table. 
The widths of the voids are measured in the direction shown by the colored lines. 
Each color represents a separately-detected void, (b) The distribution of measured 
void widths at <f> rs 0.8. Dashed line is a fit to an exponential distribution. 



In particular, it remains unclear how to correctly incorporate the partial and/or 
global pressure P. Contrary to the assumptions of the model, Smart and Ottino 
[29] used simulations to observe that over a large parameter regime, the pressure on 
small and large particles is in fact close to equal. Furthermore, Golick and Daniels 
[3"0] observed that the rate at which shear-segregation occurs (related to S) is a non- 
monotonic function of the particle size ratio, with dissimilar sizes experiencing the most 
sensitivity to the global pressure. Finally, in this Focus Issue, Fan and Hill [3T] have 
had success with a model similar to that of Gray and Thornton, with the addition 
of terms which account for gradients in kinetic stresses, as measured through velocity 
fluctuations. 

Modern experimental techniques, as well as discrete element simulations, now allow 
us to measure both the local void distribution and the trajectory of each particle. 
This provides all of the local dynamics needed to directly observe particles falling into 
voids. Fundamentally, the assumptions behind the kinetic sieving mechanisms are both 
empirically testable ("What is the difference in pressure on large vs. small particles?") 
and related to recent theories about the statistical mechanics of static granular materials. 
Using theories which describe the distribution of free space [31 @] it should be possible 
to calculate transition probabilities for either large or small grains within a bath of a 
distribution of known void sizes. Furthermore, theories which include the stresses on 
the grains [32], [331 El] would permit the direct calculation of pressure-imbalances. 

Interestingly, ensemble-approaches within the granular physics community have 
focused on the local distributions, whereby all un-occupied space is partitioned up 
and associated with the nearest particle. Two popular measures of are the radical 
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Voronoi tessellation and the navigation map [36] . However, other scientific communities 
[37] 138] 139] 140] have focused on the void-size (e) distribution, an example of which is 
shown in Fig. [3^, for a two-dimensional granular system. The void distribution is of 
primary importance for segregation problems, as a means to construct an empirical or 
analytical expression of the maximum-entropy ideas first laid out by Savage and Lun 
[22] . In Fig. ^p, the approximately exponential shape is similar to what is predicted by 
Eq. [T] The light/dark shaded regions represent the probability of finding a void space 
which is large enough for small/large particles to move into. The ratio of these two 
probabilitiess, even for the simlarly-sized particles shown here, is more than 5. Direct 
measurements of function of the local environment (£, P, <p, e) will allow us to 

determine the most relevant dependencies. 

3. Axial segregation and coarsening in a rotating drum 

In many industrially- or geophysically- relevant systems, the flow of particles is less rapid 
and collisional than those described in the previous section, and particles are instead 
densely-packed and in sustained contact with each other. Rather than strong gradients 
driving the flow, particles slowly roll and slide past each other. The paradigmatic case for 
segregation in this dense-flow regime is the rotating drum: a partially-filled cylindrical 
container containing two sizes of grains which rotates slowly around its horizontally- 
oriented axis [4T1 142] . 

When the rotating drum is tube-shaped (longer than it is thick), three sequential 
segregation stages can be distinguished. Within the first few rotations, radial segregation 
occurs: the small particles migrate towards the central axis of the drum and form a 
channel extending the whole length of the cylinder. This channel is surrounded by a 
nearly-pure layer of large particles; it forms via a kinetic sieving similar to that discussed 
in the previous section. 

During the next ten to hundred rotations, a process called axial segregation happens: 
the initial core of small particles becomes unstable to longitudinal modulations and 
finally alternating bands of large and small particles appear along the axis of the drum. 
In addition, the bands of small particles remain connected by the radially-segregated 
core of small particles as shown in Fig. [4j see also the accompanying Movie 1. 

The basic mechanism behind axial segregation is still an open question. As most 
of the grain motion happens at the free surface, it was initially thought that axial 
segregation is due to differences in the dynamic angle of repose of the two particle sizes 
[43J. Indeed there are clear differences in the surface flow profiles of the two sizes [33]. 
This idea has subsequently been expanded to a full continuum description by adding 
the relative concentration of particles [45, 46J. 

With the advent of three-dimensional imaging measurements (primarily MRI) 
it became evident that subsurface flow (or, more precisely, the migration of small 
particles in the radial core) is the mechanism which creates the axial segregation pattern 
[47] 145] 1451 loD] I5T] . This result has also been supported by both numerics (44] [52] and 




Figure 4. Axial segregation pattern in a rotating drum experiment. The tube was 
half filled with equal volumes of particles with diameters 1 mm and 0.425 mm, and 
then rotated at 15 rpm until the axial stripes had formed. Image (a) shows the outer 
appearance of the bed, while (b) highlights the connecting core between the same two 
stripes, by removing the large particles from the image. Both images are rendered with 
POV-Ray using experimental particle positions determined by X-ray tomography. A 
more complete visualization is provided by the accompanying Movie 1. 



experiments showing that mixtures of same size spheres but different angles of repose 
do not show any axial or radial segregation [53]. However, the question of the driving 
force behind the flow through the core is still not answered. While it has been shown 
that frictional interactions at the end walls can support the formation of axial bands 
|54j . axial segregation appears simultaneously everywhere in long drums, rather than at 
first in the vicinity of the end walls [551 EU [50]. Below, we suggest a new origin for this 
radial core flow. 

Finally, on timescales of hundreds to many thousands of rotations, these axially- 
segregated bands undergo coarsening. During this stage (see Fig. |5^i), the large-particle 
bands grow together and merge, while the number of small-particle bands continually 
decreases. Finger et al. [56] used carefully-prepared samples containing only two bands 
of small particles to show that it is typically the band of smaller volume which vanishes. 
The process by which this happens corresponds to a continuous flow of particles through 
the connecting radial channel, rather than a diffusive random walk [571 HU EHJ [52j [51] . 
The only known theoretical approach capable of explaining the coarsening dynamics is 
the continuum theory laid out in [4"5l H6] . but it ignores the dynamics inside the core. 

In this review, we suggest that a unified explanation may be able to account for 
both the axial segregation and the subsequent coarsening. It is helpful to first consider 
the often- used analogy of a granular fluid [SB]: axial segregation is morphologically 
similar to the Rayleigh-Plateau instability in ordinary fluids, whereby a cylindrical fluid 
column breaks up into individual drops. The instability is driven by the reduction in 
surface energy provided by the smaller surface area of the droplet state. 

In the case of granular fluids, such an explanation would depend on the existence of 
a effective surface tension a between the two fluids (small-particle and large-particle) ||] 

J Such a surface tension is distinct from the effective surface tension which arises due to small cohesive 
forces, recently found in freely falling granular streams [551 15U] . 
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Figure 5. Coarsening dynamics in a rotating drum experiment, (a) Space-time plot 
of the evolution inside a long rotating drum. Dark areas correspond to bands of small 
glass spheres (0.55 mm diameter), bright areas consist of large glass spheres (1.5 mm). 
The drum rotates at 15 rpm. Reprint from Finger et al. [55] • (b, c, d) Vertical cross 
section through the middle of a rotating drum initially containing two bands of small 
particles. Images are taken with x-ray tomography after 100, 9500, and 24300 rotations 
with 20 rpm. Courtesy Tilo Finger. 



Such a surface tension, to be described below, would also provide an intuitive explanation 
for the coarsening dynamics: it is a consequence of the system continuing to reduce its 
interfacial area by merging several small-particle bands into a single, extended band. It 
is worth noting that the kinetic sieving mechanism discussed in the last section provides 
a stabilizing mechanism at the necks which connect the individual bands/droplets. Only 
because the necks don't rupture as in an ordinary fluid (see Fig. [5|Dcd), the small-particle 
bands can continue to exchange "fluid" till they have merged into one extended stripe. 

How might such an interfacial free energy arise between two granular fluids? 
In dense granular system, the interparticle interactions take the form of Hertzian 
normal forces and frictional tangential forces. Neither of these forces depends strongly 
on particle size. Therefore, the free energy will mostly be determined by entropic 
contributions. The enduring or sliding contacts make the system quite different from 
a classical hard sphere systems, where the configurations with even only two particles 
touching are of measure zero. Therefore we will turn to the Edwards ensemble [5j |6] for 
a framework in which to discuss the relevant configurational entropies. 

Edwards and coworkers suggested the utility of a configurational entropy 

Sconf = A In Q, (3) 

where Q is the number of all possible mechanically stable configurations under 
constant volume fraction and boundary stresses, and A the as-yet-undetermined granular 
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equivalent to the Boltzmann constant. While there are a number of open questions, for 
instance the role of ergodicity, a growing community has started to work on such a 
statistical mechanics approach based on the Edwards ensemble [61J. 

For the purposes of determining an interfacial tension, the salient observation is 
that the number of possible packings of monodisperse spheres does not depend on their 
diameter at given volume fraction. Therefore, S con f/N (entropy per particle) is also 
independent of the diameter for a pure granular fluid. However, at an interface between 
two sizes of particles, the number of mechanically stable configurations can be expected 
to depend on the diameter ratio. As the ratio between the two sizes deviates from 
unity, there will be extra void space at the interface which can not be effectively used 
to create packings. This is similar to granular interfaces with flat walls [62| 163] . This 
unusable void space will result in a reduced value for S con f, as compared to the value 
for a single-component fluid, and will therefore provide an effective surface tension. 

Such an effective surface tension would also explain why the narrower of two 
small-particles bands is typically the one which vanishes, as shown in Fig. [5] The 
narrower band has a less favorable surface-to- volume ratio, resulting in fewer accessible 
configurations as compared to the wider band. Therefore, particles within the connecting 
channel will find themselves more often in configurations further away from the center 
of the narrower band. This results in a net flux from the narrower to wider band. 

This analogy to interfacial tension remains, at present, a plausibility argument 
rather than a proof of the existence of an effective surface tension. We see two potential 
ways of testing this hypothesis. First, one could manually prepare a radially-segregated 
system in which the initial cylindrical interface between the two fluids is as smooth as 
possible. A series of X-ray tomographs covering the first few rotations will permit 
measurements of the evolution of the interface. If the evolution is indeed driven 
by an effective surface energy, a Fourier-decomposition of the interface will reveal a 
characteristic range of exponentially growing modes. In the Rayleigh-Plateau scenario 
[6%] the wavelength of the fastest growing mode is independent of the value of a . 
However, the growth rate should depend on a, which in turn can be modified by changing 
the size ratio of the particles. 

Second, one could use capillary wave theory to test for the existence of an effective 
surface tension [651 [66J . This would require a sheared setup in which the large particles 
initially form a layer on top of the small particles. As this sequence of layers is 
stable to kinetic sieving, the externally-imposed shear will create small fluctuations 
of the interface, which can again be captured by x-ray tomography. The mean square 
amplitude (A 2 ) of the Fourier modes of the interface will depend on the wavenumber 
k and the interface sidelength L. In the presence of an effective surface tension, (A 2 ) 
should scale like 

^ ~ z^f • < 4) 

which can readily be tested. 



REFERENCES 



11 



4. Conclusion 

In this review, we have described segregation in two specific granular systems of 
medium to high density. In both cases, an approach based on entropic or probabilistic 
arguments seems promising explanations of observed phenomena. We believe that 
this likely holds true for most such granular systems: the kinetic energy does not 
influence the dynamics strongly, nor is it likely that there are significant differences 
in contact forces between different particle-species. Instead, volume exclusion and 
geometrical constraints have a strong, yet mostly unexplored, influence on the number 
of available configurations (entropy). Even though the myriad ways of external driving 
such systems are presumably far from ergodic, it is the accessible phase-space volume 
which determines the fate of the system. 

It is interesting to note that while the question of segregation was pivotal to Edwards 
and coworkers when suggesting their statistical mechanics approach (5J [6] , it has only 
recently begun to see application. The work to date has primarily concerned simulated 
granular systems [67J EU ESI [70] , where particle-scale measurements of velocities, free 
volumes, and forces have been possible. With the increasing availability of experimental 
methods such as photoelastic force measurements, MRI, and X-ray tomography, the time 
is ripe to move forward understanding the particle-scale origins of these eye-catching 
bulk phenomena. 
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